Microscopic theory of surface-enhanced Raman scattering in noble-metal 

nanoparticles 
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We present a microscopic model for surface-enhanced Raman scattering (SERS) from molecules 
adsorbed on small noble-metal nanoparticles. In the absence of direct overlap of molecular or- 
bitals and electronic states in the metal, the main enhancement source is the strong electric field 
of the surface plasmon resonance in a nanoparticle acting on a molecule near the surface. In small 
particles, the electromagnetic enhancement is strongly modified by quantum-size effects. We show 
that, in nanometer-sized particles, SERS magnitude is determined by a competition between several 
quantum-size effects such as the Landau damping of surface plasmon resonance and reduced screen- 
ing near the nanoparticle surface. Using time-dependent local density approximation, we calculate 
spatial distribution of local fields near the surface and enhancement factor for different nanoparticles 



I. INTRODUCTION 

Surface-enhanced Raman scattering (SERS) has been 
one of the highlights of optical spectroscopy in metal 
nanostructures during past 25 yearsi Recent inter- 
est in SERS stems from the discovery of extremely 
strong single-molecule SERS in silver nanoparticle 
aggregates)^ as well as from nanoparticle-based applica- 
tions such as, e.g., biosensors- that rely on sensitivity of 
SERS to small concentrations of target molecules. The 
main mechanism of SERS has long been known as elec- 
tromagnetic (EM) enhancementi^iSil of dipole moment 
of a molecule by the strong local field of surface plas- 
mon (SP) resonance in a nanoparticle. EM mechanism 
is especially effective when a cluster of nanoparticles is 
concentrated in a small region ("hot spot" ) .guLifiiil A 
combined effect of SP local fields from different particles 
acting on a molecule trapped in a gap can result in a 
giant (up to 10 14 ) enhancement of the Raman scattering 
iilfiilfiill Other mechanisms contribut- 
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ing to SERS can involve electron tunneling between a 
molecule and a nanoparticle^* 

The conventional description of EM enhancement is 
based on classical Mie scattering theory&£ The dipole 
moment of a molecule at distance ro form a particle cen- 
ter is enhanced by a factor ~ a p (ui)/rQ, where a p = 
i? 3 [(e — l)/(e+ 2)] is the particle polarizability, R is its 
radius, and e(u>) is metal dielectric function. The far-field 
of molecular dipole, radiating at Stokes-shifted frequency 
u) a , is, in turn, comprised of direct and Mie-scattered 
fields. The latter contributes another factor ~ a p (ui s )/i"Q, 
so that the total field enhancement is ~ a p (u>)a p (oj s )/rQ 
and Raman crossection is proportional to |a p | 4 /rQ 2 . At 
frequencies close to the SP pole in ev p , this enhancement 
can reach ~ I0 6 . Note that, within classical description, 
the dependence of SERS on nanoparticle size, coming 
from geometrical factor in a, is weak if the molecule is 
sufficiently close to nanoparticle surface. 

The classical approach is valid for relatively large 



nanoparticles, where the effect of confining potential on 
electronic states is negligible. For nanoparticle sizes 
R < 5 nm, the lifetime of SP is reduced due to the Lan- 
dau damping by single-particle excitations accompanied 
by momentum transfer to the surface. 1 ^ This results in 
a broadening of SP resonance peak by the amount of 
level spacing at the Fermi energy, 7 S ~ vp/R (vf is the 
Fermi velocity) , and in the corresponding decrease of the 
SP field amplitude. For not very small nanoparticles, 
the effect of Landau damping on SERS can be treated 
semiclassically* by incorporating the quantum-size cor- 
rection 7 S in the Drude dielectric function of metal. Re- 
cent resonance fluorescence measurements on small gold 
particlesj*2*2ii2* however, indicate strong deviations of 
the plasmon-induced enhancement from that predicted 
by semiclassical models^ For clusters with electron num- 
ber N < 100, SP lifetime is reduced to several fs and 
SERS is diminished. Note that some enhancement of 
Raman signal due to single-particle resonances remains 
even in small clusters containing several atomsi^* 

In this paper, we study SERS for small nanoparticles 
several nanometers in diameter, i.e., in the intermedi- 
ate regime between classical particles and small clus- 
ters. Our chief observation is that, in this crossover 
regime, SERS magnitude is determined by interplay be- 
tween several competing quantum-size effects, including 
the aforementioned SP Landau damping as well as the 
modification of electron screening in the surface region. 
The latter produces an opposite trend towards a rela- 
tive increase of SERS when the molecule is located in 
a close proximity to the metal surface. The underlying 
mechanism is related to different effects that the confin- 
ing potential has on rf-band and sp-band electron states. 
Namely, the deviation of potential well from the rectan- 
gular shape gives rise to a larger effective radius for the 
higher-energy sp-electrons^* This effect is further am- 
plified by the spillover of sp-band electron density, due 
to tunneling into the barrier, as contrasted to the es- 
sentially step-like density profile of localized d-electrons. 
As a result, in a surface layer of thickness I — 2 A, the 
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G?-electron population is diminished and, hence, the in- 
terband (i.e. due to d-sp transitions) polarizability is 
strongly reducedi 2 ^ 2 ^ The resulting reduction of screen- 
ing in the surface region leads to a greater strength of 
electron-electron interactions in the sp-band that was ob- 
served, e.g., in a faster, as compared to bulk metal, elec- 
tron relaxation in Ag nanoparticles measured using ultra- 
fast pump-probe spectroscopy^ It is, therefore, natural 
to expect that such underscreening should lead to addi- 
tional enhancement of local field outside of the nanopar- 
ticle. Note, however, that here the effect of spillover is 
twofold: while it leads to an increase of volume fraction 
of underscreened region and hence to stronger local field, 
especially in small nanoparticles, it can also, by itself, 
have an opposite effect by smearing out the otherwise 
sharp classical boundary. 28 Thus, in small nanoparticles, 
the enhancement magnitude is determined by a delicate 
interplay of competing quantum-size effects, and must, 
therefore, be described within a consistent microscopic 
approach. Such an approach, based on time-dependent 
local density approximation (TDLDA) ^ is developed in 
this paper. 

The paper is organized as follows. In Section^] we de- 
rive the general expression for polarizability of molecule- 
nanoparticle system. In Section IIIII we derive a self- 
consistent system for local potential that determines the 
enhancement factor. In Section Hvl we present the results 
of our numerical calculations. Section [V] concludes the 
paper. 



II. POLARIZABILITY OF 
MOLECULE-NANOPARTICLE SYSTEM 

We start with formulating SERS in terms of quantum 
transitions in the interacting molecule-nanoparticle sys- 
tem. We assume that the molecule is located at distance 
ro from the nanoparticle center and that the overall sys- 
tem size is much smaller than radiation wavelength so the 
retardation effects can be ignored^ In the absence of di- 
rect electron tunneling^ the interactions within excited 
molecule-nanoparticle system are caused by nonradiative 
transitions accompanied by energy transfer between a 
molecule and a nanoparticle, similar to Forster transfer 
in two-molecule systems^ Namely, an electron-hole pair 
can nonradiatively recombine by transferring its energy 
to SP (and vice versa) via dynamically-screened Coulomb 
interaction^ Feynman diagrams of processes contribut- 
ing to polarizability of molecule-nanoparticle system, a, 
are shown in Fig. ^ (a) incident photon with energy u> 
is absorbed by the molecule and reemitted with Stokes- 
shifted energy u> s ; (b) after absorbing a photon, excited 
molecule nonradiatively recombines transferring its en- 
ergy to SP in the nanoparticle, which emits a photon; 
(c) SP, excited by incident light, transfers its energy to 
the molecule, which emits a photon; and (d) after energy 
transfer from SP to molecule, the latter transfers the en- 
ergy back to SP, which emits a photon. Correspondingly, 



(a) 

(a) 



ox 



(b) 



CO 



CO. 



(c) (d) 



FIG. 1: Nonradiative processes contributing to Raman scat- 
tering from a molecule-nanoparticle system. 



the system polarizability is (in operator form) 
a = a + aUII + UUa + UUallU, 



(1) 



where a is molecular polarizability, II is density-density 
response function of a nanoparticle in medium, and U 
is the Coulomb potential. The Raman polarizability is 
obtained by calculating the matrix element of a between 
incoming and outgoing photon states with energies u) and 
lu s , respectively. The interaction of the molecule with the 
nanoparticle involves matrix elements {e\(f>(r)\g} , where 
\g) and |e) stand for the molecule ground and excited 
electronic bands, respectively, and 

0(w,r)= /"drKfra^fr-rOnfw.ri.raJ^oCra) (2) 

is the nanoparticle response to external photon potential, 
4>o( r )- Since the length scale of 0(r) is much larger than 
the molecule size, we have (e\<j>(r)\g) w /xV(/>(i"o), where 
IX is the dipole matrix element of corresponding molecu- 
lar transition^ The averaging over random orientations 
of fj, can be accounted by assuming isotropic Raman po- 
larizability tensor a. The nanoparticle contribution then 
factors out, a = otM, with 



M = 1 + — * 
El 



E • V<£(w, r ) + E • VcKw s , r ) 

+V^,r ).V^ s ,r ) , (3) 



where Eo is the electric field in the absence of nanopar- 
ticle. For incident field, Ej, polarized along the z-axis, 
<pa = —eEiT cos 0, we have Eo = Ej/e TO , where e m is the 
dielectric constant of medium. 



III. ENHANCEMENT FACTOR 

To evaluate the local potential <p(u},r) within TDLDA 
approach, we present it in the form 



, I l3 , Sn(u>, r') 
(^,r)=e 2 / dV v _' r > 



(4) 



where the induced density, 

5n(r) = J dr'Il(r,r')(t>o(r') = 5n s (r) + 5nd(r) + 5n m (r), 

(5) 
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contains contributions from sp-electrons, d-electrons, and 
surrounding medium, respectively (hereafter we suppress 
frequency dependence). 

We adopt the two-region model that combines a 
quantum-mechanical description for sp-band electrons 
and phenomenological treatment <i-electrons with bulk- 
like ground-state density n d in the region confined by 
R d < R. 2 ^ This model has been used for calcula- 
tions of polarizabilities of small Ag nanoparticles and 
clusterS(2a2a but it remains reliable for relatively large 
electron numbers, N > 1000. The induced density of 
sp-band electrons is determined from TDLDA equation 



with a = Rd/R- The kernel K (r, r'), relating the induced 
potential and density of sp-electrons, is given by 

K (r, r ) = u(r, r') 

~f3(r/R d ) u(R d ,r') -2aX m u(R,r') \d/r) 

+(3(r/R) [u{R, r') - a 2 X d u{R d , r')] X m /n, (12) 
where u(r, r') — 4irr < /3r 2 > is the dipole term of Coulomb 



potential expansion and (3{x) = x 



• 1)- 2x0(1 -ac). 



With decomposition ©, the TDLDA equation © takes 
the form 

1 



Sn s (r) = / dVi>(r,r') <S>{r') + V^[n{r')]5n s {r') , (6) Sn s (r) = / dr'r' 2 P s (r, r')-f- fatf) + 6<f> d (r') 



e(r') 



where $ = 0o + 4> is the full potential, P s (r, r') is 
the polarization operator for noninteracting sp-electrons, 
V^[n(r')] is the (functional) derivative of the exchange- 
correlation potential and n(r) is the ground-state elec- 
tron density. The latter is obtained in a standard way by 
solving Kohn-Sham equations. To close the system, we 
need to express the full potential <£>(r) via Sn s (r). This 
is accomplished be relating Sn d (r) and Sn m (r) back to 
$(r) as 



e 2 Sn d (r) = V[ W (r)V$(r)], 
e 2 6n m {r) = V[x™(r)V$(r)] , 



(7) 



where Xd( r ) = [( e d — l)/4vr] 0(Rd — r) is the inter- 
band susceptibility with the step function enforcing the 
boundary conditions and, correspondingly, Xm( r ) — 
\(e m — l)/47r] 6(r — R) is the susceptibility of surrounding 
medium. The derivation is given in Appendix. The final 
result is conveniently expressed in terms of expansion of 
$(r) and 8n{r) in spherical harmonics. For non-resonant 
Raman scattering, only dipole (L = 1) terms contribute 
to the local field. The final expressions for dipole compo- 
nent <!>(r) can be presented as a decomposition (see the 
Appendix) 



where 4>o(r) 



e(r) 



5cj) d (r) + 5<f> s {r)], 



(8) 



-eEi.r, 



b d (r) = (3(r/R)MR) A m (l - a 3 X d )/n 

-P(r/R d )MRd) A d (l - 2X m )/n (9) 



and 



(r) = / dr'r' 2 K(r,r')5n s (r'). 



(10) 



Here e(r) — (e d , 1, e m ) for r in the intervals [(0, R), 
(Rd,R), (R,oo)], respectively, and 



X,j 



e d + 2 : 



X, 



2e r , 



V =l-2a 3 X d X m , (11) 



dr"r" 2 K(r',r")6n s (r") 



+V^r')Sn s (r') 



(13) 



Note that <&(r) is continuous at r — R d , R. 

Equations ffSf) — (| 1 3|) determine self-consistently the spa- 
tial distribution of local potential near small noble-metal 
nanoparticles. Here 5<j) d {r) is the induced potential due 
to (^-electrons and surrounding medium. Their effect on 
the sp-electron potential, 5<j) s (r) , is encoded in the kernel 
K(r, r 1 ). For e d = e m = 1, we have K(r, r') = u(r, r') and 
54>d(r) — 0, recovering the case of simple metal particles 
in vacuum. 

If the molecule is not too close to the surface (d > 
lA), i.e., there is no significant overlap between molecular 
orbitals and electronic states, then Eqs. Q and ijSJ 113|) 
yield 

1 eE- 
<5$(r ) ee — - [S(f> d (ro) + ^M] = V a v > ( 14 ) 

where a p {u>) is the nanoparticle polarizability. The ex- 
pression for a p is given in Appendix. Field enhancement 
coefficient M, Eq. Q, then takes the form 



= 1 + (1 + cos 2 fl ) 



M = 



+(l + 3cos 2 0o) ^#M , (15 ) 

In this case, enhancement retains the same functional de- 
pendence on particle polarizability as in classical theory^ 
however, a p is now determined microscopically. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Below we present our results for SERS enhancement 
factor \M\ 2 for Ag nanoparticles in a medium with di- 
electric constant e m = 1.5. Calculations were carried for 
number of electrons ranging from N — 92 to TV = 3028, 
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corresponding to particle diameters in the range D m 
1.4 — 4.5 nm. To ensure spherical symmetry, only closed- 
shell "magic numbers" were usedi£& For such sizes, the Ag 
band-structure remains intact. The ground state energy 
spectrum and wave-functions were obtained by solving 
the Kohn-Sham equations for jellium models with the 
Gunnarsson-Lundqvist exchange-correlation potential; 37 
the interaction strength was appropriately modified to 
account for static rf-band screening. The ground state 
density of sp-band electrons, n(r), exhibits characteristic 
Friedel oscillations, while the spatial extent of spillover 
is ~ 2 a.u. (see Fig. [2J - 

These results were used as input in the numerical so- 
lution of TDLDA system The molecule was 
located along the z-axis (9o — 0) at a distance d — 5 a.u. 
from effective boundary with radius R — r s N x ^ , where 
r s = (47rn/3) -1 / 3 (r s — 3.0 a.u. for Ag), ensuring no 
direct overlap with the nanoparticle. In the calculation 
of optical response, the experimental data for e^iuS) in 
Ag was used 3 ^. We also assumed that the Stokes shift is 
much smaller than 7 S ~ vp/R (j s ~ 0.5 eV for D = 3.0 
nm) and ignored the difference between ui and ui s . 

The calculated absorption spectra for different 
nanoparticle sizes are shown in Fig. [3] For medium di- 
electric constant e m = 1.5, the position of SP resonance 
at ~ 3.2 eV is well below of the interband transition onset 
at ~ 4.0 eV. As expected, with increasing size the peak 
width is reduced due to a weaker Landau damping of SP. 
In order to illustrate the role of interband screening, we 
also show the results of calculations with A = 0. Note 
that decreasing the surface layer thickness by 1.0 a.u. 
only somewhat reduces the underscreening; the main ef- 
fect still comes from the larger spatial extent (about 2 
a.u.) of sp-band electron spillover. The effect of un- 
derscreened surface region on the absorption is two-fold. 
The large contrast ratio of e c i in the bulk and surface re- 
gions [€d(uJsp) ~ 5] leads to a lower, as compared to bulk, 
average value of interband dielectric function in the rel- 
evant frequency region. As a result, the peak position 
for A = 1.0 a.u. is slightly blueshifted with respect to 
that for A = (Sui sp ~ 0.05 eV). At the same time, for 
A = 1.0 a.u., the peak amplitude is larger although the 
resonance width stays unchanged. The stronger absorp- 
tion for A = 1.0 a.u. is caused by a weaker screening of 
the SP electric field, that determines the peak oscillator 
strength, in the surface region. 

The calculated local field, E, at resonance frequency 
is plotted in Fig. 0] vs. molecule-nanoparticle distance, 
d = ro — R. The gradual rise of field magnitude on 
the length scale of electron spillover replaces the discon- 
tinuity (for a, e m ^ 1) of classical field at the sharp 
boundary. It can be seen that while at A = 1.0 a.u. 
the field amplitude is larger than for A = 0, the differ- 
ence decreases for larger d. Correspondingly, the effect 
of interband screening on Raman signal enhancement, 
\M\ 2 oc |-E| 4 , is substantial only when the molecule is 
located close to the nanoparticle. 

In Fig. \5\ we plot the enhancement factor \M\ 2 as a 




r(a.u.) 

FIG. 2: Calculated ground state density n(r) for Ag nanopar- 
ticles is shown for closed-shell electron numbers N — 92 (a), 
N = 556 (b), and N = 3028 (c). 

function of incident light energy for different nanopar- 
ticle sizes. Note that the asymmetric shape of the en- 
hancement peak, as compared to the absorption peak, is 
because the former is determined by the absolute value of 
polarization, rather than its imaginary part. The general 
tendency is a decrease of SERS for smaller nanoparticles 
with the enhancement factor reaching only \M\ 2 ~ 100 
for the smallest nanoparticle size, D ~ 1.4 nm. This de- 
crease is mainly related to the Landau damping of SP: 
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FIG. 3: Calculated absorption spectra for Ag nanoparticles 
with A = (dashed line) and A = 1.0 a.u. (solid line) is 
shown for electron numbers N = 92 (a), N = 556 (b), and 
TV = 3028 (c). 



at resonance energy, we have E rx I/7 ~ R/vp yield- 
ing \M\ 2 oc R 4 . The comparison of results for A = 
A = 1.0 a.u. shown in Figs.JSfa) and|S{b), respectively, 
shows that reducing the thickness of surface layer even 
by 1.0 a.u. substantially affects the enhancement. The 
role of interband screening is most visible in the depen- 
dence of SERS on nanoparticle size plotted in Fig. 
at resonance frequency for nanoparticle diameters in the 
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FIG. 4: Calculated local field for several nanoparticles sizes 
vs. distance from classical boundary for A = and A = 1.0 
a.u. 
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FIG. 5: Calculated enhancement factor \M\ 2 for different 
nanoparticle sizes is shown with A = (a) and A = 1.0 a.u. 
(b). 



range 1.4-4.5 nm. With decreasing size, the enhance- 
ment factor drops by an order of magnitude, while over- 
all enhancement is reduced for A = 0. Importantly, the 
latter effect is more pronounced for smaller nanoparti- 
cles; for A = 1.0 a.u., the enhancement factor is larger 
by 25% for N = 3028 but by 35% for N = 92 than for 
A = indicating a more important role of screening in 
smaller nanoparticles due to larger surface-to- volume ra- 
tio. Thus, the proper account of screening gives a slower 
(as compared to semiclassical models) decrease of the en- 
hancement as nanoparticles become smaller. Note, fi- 
nally, that by keeping A constant for different nanopar- 
ticle sizes, we somewhat underestimated the screening 
contribution to SERS; indeed, for smaller nanoparticles, 
underscreened layer is thicker due to stronger deviations 
of the confining potential from rectangular shape. 
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frequency vs. nanoparticle size is shown with A = and 
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CONCLUSIONS 



We developed a microscopic model for surface- 
enhanced Raman scattering in noble-metal nanoparticlcs. 
Our approach incorporates, in a unified manner, all rele- 
vant quantum-size effects that determine the magnitude 
of Raman signal enhancement for nanometer-sized par- 
ticles. While the Landau damping of surface plasmons 
leads to a general decrease of the enhancement for small 
particles, this trend is partially offset by the reduction of 
interband screening in the metal boundary region. The 
additional enhancement of local field is substantial only 
in a close proximity to the metal surface, and it is more 
pronounced for smaller nanoparticles where the electron 
density profile deviates strongly from the classical shape 

As a final remark, we considered here a nonresonant 
Raman scattering, i.e. the case when the excitation en- 
ergy of a molecule is much larger than the SP energy 
uj sp . In this case, due to small single- molecule polariz- 
ability, only linear (in molecule) response needs to be 
considered and SERS is dominated by a single back and 
forth energy transfer within molecule-nanoparticle sys- 
tem [see Fig.^d)]. In the resonant case, multiple energy 
transfer processes give rise to a nonradiative width of 
the molecular levels and, in general, to a reduction of the 
enhancement^ For small molecule-surface separations, 
these nonradiative processes are dominated by surface- 
enhanced electron-hole pair generation in the metal. 
Even in simple metals, these processes are enhanced due 
to a reduced intraband (Thomas-Fermi) screening near 
the boundary4Mi The issue of nonradiative decay for 
noble-metal nanoparticles, where interband screening ef- 
fects become important, will be addressed in a future 
publication. 



Acknowledgments 

This work was supported by NSF under Grants No. 
DMR-0305557 and NUE-0407108, by NIH under Grant 
No. 5 S06 GM008047-31, and by ARL under Grant No. 
DAAD19-01-2-0014. 



APPENDIX: DERIVATION OF LOCAL 
POTENTIAL 

The full self-consistent potential <&(]") = <j>(r) + (f>o(r) 
sutisfi.es Poisson equation 

$( W ,r)=0 o (r) + e 2 |dV^^, (A.l) 

where the induced density si comprised of sp-band, ri- 
band and medium contriburions, <5rt(r) = Sn s (r) + 
Snd(r) -(- 5n m (r). Using Eq. JJJ for Sn m and Snd and 
integrating by parts, Eq. (JA.lfl takes the form 



e(r)$(r) = fo(r) + e 2 / d 6 r 



3 ,Sn s (r') 



4tt 



f dS/v'—l . V6{R d - r)cfr(r') 
J |r — r'| 



+ f d z r >sj'—±— . v'0(r - i?)$(r'), (A.2) 

An J |r — r'| 

where e(r) — (e d , 1, e m ) for r in the intervals [(0, R), 
(Rd, -R), (R, oo)], respectively. Since the source term has 
the form 4>o(r) = 4>o(r)cos9 = —eE\r cos 8, we expand 
$ and 5n s in terms of spherical harmonics and, keeping 
only the dipole term (L = 1), obtain 



e(r)$(r) = 4>o(r) + e 2 J dr'r' 2 u(r, r')5u s (r') 
e d - 1 2 9u ( r > R d) 



An 



Air ~ Rd ' dR d 

e m - 1 d2 du(r, R) 
-Jx 



dR 



HR d ) 
*(R), (A.3) 



where 



47T 

u(r, r) = — 
v ' ' 3 



^0(r-r') + -^ e(r'-r) 



(A.4) 



is the dipole term of the radial component of the Coulomb 
potential. The above equation can be simplified to 

e(r)$(r) = 0(r) - ^±(3(r/R d )$(R d ) 

+^-^P(r/R)$(R), (A.5) 



where (3(r/R) = (3i? 2 /47r) [du(r, R)/dR] is given by 

P(x) = x- 2 6(x - 1) - 2x0(1 - x), (A.6) 
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and we introduced a shorthand notation 

4>{r) = <p (r) + e 2 J dr'r' 2 u(r, r')Sn s (r r ). (A.7) 

The boundary values of can be obtained by matching 
$(r) at r = R d , R, yielding 

(e d + 2)${R d ) + 2a(e m - 1)§(R) = 3$(R d ), 

(e d - l)aH(R d ) + (2e m + 1)$(JJ) = ty(R), (A.8) 

where a = R g/R. Substituting <&(R d ) and $(i?) back 
into Eq. I|A.5J) . we arrive at 

e(r)$(r) = 0(r) - f3(r/R d ) ^ U{R d ) - 2a\ m 0(R) 

f] L 

+0(r/R) — \$(R) - a 2 \ d 0(R d )] , (A.9) 
rj L J 

where the cofhcients A are given by Eq. (|llf) . Separating 
out 6n s -dependent contribution, we arrive at Eq. (JSJ. 

The expression for nanoparticle polarizability a p (u>), 
can be obtained from the large-r asymptotics of induced 
potential, Eq. JQ. From Eqs. we find for 



r > i? 



<50 d (r) = — ^ a d , (50 s (r) = — ^ a s , 



(A.10) 



with 



3 (1 — A m )a 3 A ( j — A r 



1 - 2a 3 A d A r , 



4-7T 



drV /3 5n s (r') 



A d - A m + A d (l - 2a 3 ) ^, ,j. 



1 - 2a 3 A d A m jo 
(1 — A m )a 3 Ad — A m 3 



aYV' 3 on s (r') 



1 - 2a 3 A d A m 



i? 3 / dr'Sn s (r' 

R 



m JR 



.(A.ll) 



The nanoparticle polarizability is a p = a d + a s . Note 
that for frequencies below interband absorption onset, 
the d-band/medium contribution a d is real. 
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